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1 Introduction 

In the last years many experimental and theoretical efforts have been devoted to search Quark 
Gluon Plasma (QGP) and/or collective effects^] in heavy ion collisions |l], ^, ||, ||, |5], |6|, as a 
tool to study the nonperturbative aspects of Quantum Chromodynamics (QCD). To achieve 
this goal several signatures of the phase transition(s) from confined to deconfined quarks and 
gluons and of chiral symmetry restoration have been proposed. 

The finding of three of the proposed [[?], ||, ^] signals at the Super Proton Synchrotron 
(SPS) at CERN has originated great excitement and debate in the scientific community in 



this field. In fact, an abnormal suppression of J/tf) [10] and a strong enhancement of strange 



baryons and antibaryons [11] have been observed in central^ Pb-Pb collisions, compared with 



those measured in collisions between lighter projectiles and targets. Also an enhancement in 
the dilepton spectrum for dilepton masses below 0.8 GeV/c 2 has been seen in Pb-Au collisions 



[12]. Whether or not these three experimental observations are really an unambiguous proof 



of the existence of a QGP is still an open question gj| |l|, [T|, ||, 0, H, H, [p H, due 



to their possible explanation using more conventional, but still interesting, physics. In any 
case, it is expected that the forthcoming heavy ion experiments in the Relativistic Heavy 
Ion Collider (RHIC) at BNL and the Large Hadron Collider (LHC) at CERN will clarify 
definitively the point| . 

In spite of the work done these last years there are many fundamental aspects of the 
physics of heavy ions at high energies which are not clear at all. Fundamental questions like, 
e.g.: 

• Is particle rapidity density proportional to the number of participant nucleons or to the 
number of elementary nucleon-nucleon collisions? 

• Which is the physical explanation of the SPS particle correlation data? 

• How large will particle multiplicities be at RHIC and LHC? 

are answered in a very different way by several models, all of them claiming to agree with 
the existing experimental data. 



1 Collective effects are those considered to explain the event, which go beyond the superposition of elemen- 
tary nucleon-nucleon collisions. 

2 By central collision we mean a head-on one, in which most of the matter of the lightest nucleus participates. 

In practice different criteria are used, both theoretically (upper bound in impact parameter, minimum number 

of participant wounded nucleons,. . .) and experimentally (percentage of the cross section, lower bound in the 

number of charged particles,. . .). 

3 RHIC and LHC will provide center of mass energies of 200 GeV and 5.5 TeV per nucleon respectively, to 

be compared with ~ 20 GeV per nucleon at the SPS. 
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Referring to the last question, in Fig. [I] it is shown the pseudorapidity distribution of 
charged particles from different models for central Pb-Pb collisions at a beam energy of 3 
TeV per nucleon; this plot has been taken from the ALICF^j Technical Proposal [22] done by 
the ALICE Event Generator Pool in December 1995. The results according to Monte Carlo 
codes of several models show large differences at central pseudorapidity. Indeed, between the 
String Fusion Model (SFM) || and the VENUS || or SHAKER ||| codes there is a factor 
larger than 4 at rj = 0, while the difference in the fragmentation regions (\rj\ > 5) is smaller. 
At RHIC energy^] the difference is about a factor 2 - most models give results in the range 
700 -T- 1500. 

These uncertainties in one of the most elementary aspects of the collision, may leave 
us uncomfortable regarding the necessity to keep under control the conventional physics 
of heavy ions to clearly distinguish the signatures of QGP and/or collective effects in the 
proposed observables. Needless to say, from the experimental point of view it is crucial to 
know whether there will be 2000 or 8000 charged particles per unit rapidity in central Pb-Pb 
collisions at the LHC for the design of the detectors. For these reasons we review in this paper 
charged particle central rapidity density predictions of different models for central collisions 
between the largest nuclei that will be available at RHIC and LHC, discussing the origin of 
the differences among the results. 

According to their origin, models can be classified into three categories. On the one hand, 
some models like Dual Parton Model (DPM) |28| , p9j ], its Monte Carlo implementation 
DPMJET Quark-Gluon String Model (QGSM) §|, FRITIOF @], SFM || H,||, 



Relativistic Quantum Molecular Dynamics (RQMD) |36|, 37, 38], Ultrarelativistic Quantum 



Molecular Dynamics (UrQMD) |3g [Ij, VENUS gj or its new version NEXUS g|, and 
LUCIAE [^] mainly pay attention to the soft part of the collision (there is no need of a hard 
perturbative part at SPS energies). The hard part in some of these models is included adding 
to the elementary soft cross section the jet one, as an input for the eikonalized cross section. 
On the contrary, other models like the Heavy-Ion Jet Interaction Generator (HIJING) 



[43], Eskola et al. [44], and Geiger and Miiller [45] are mainly focused to the hard part. They 
compute the number of minijets or partons with transverse momentum larger than a given 



p±o > 1-^2 GeV/c. These hard partons are taken [46] as the starting point of an evolution and 
expansion previous to hadronization (for discussions on this point see for example |47], |48|]), 
A soft part, extracted from the SPS data, is added with an energy dependence taken from 
some model. 



4 ALICE (A Large Ion Collider Experiment) is the approved detector at the LHC fully dedicated to Heavy 
Ion Physics. 



J Un updated set of predictions for RHIC can be found in |ir| . 
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A third kind of models are the statistical and thermodynamical ones ||. In these 
models the main predictions refer to ratios between different kind of particles and not to 
absolute values of each kind. Usually to get absolute rapidity densities the volume at freeze- 
out has to be specified. The volumes used are ~ 3600 and ~ 14400 fm 3 , giving charged 
particle densities at midrapidity of 1200 and 8000 at RHIC and LHC energies respectively^ . 

The plan of this review will be the following: After this Introduction, in the next Section 
the DPM and DPMJET Monte Carlo code will be discussed in some detail, introducing several 
concepts which will also be used in the other models. In Sections 3, 4, 5 and 6 the SFM, 
RQMD, HIJING and Perturbative QCD (PQCD) and Hydrodynamical models respectively, 
together with their predictions for charged particle densities at midrapidity, will be briefly 
reviewed. Other models will be discussed in Section 7. Afterwards, in Section 8 we will argue 
on percolation in heavy ion collisions, and in Section 9 possible implications for Cosmic Ray 
Physics will be commented. In the last Section the different results will be compared and 
some discussions presented. 



2 The Dual Parton Model and the DPMJET Monte Carlo 
code 



The DPM is a dynamical model for low p± hadronic and nuclear interactions, based on the 



large N expansion of QCD with N c /Nf fixed [49]. The dominant lowest order configuration 
in p-p scattering at high energy consists in the production of two strings between valence 
constituents, of type (qq) v — q v , see Fig. ^. 

There are also more complicated terms, corresponding to higher order diagrams in the 
large N expansion, involving 4, 6,. . . strings. These extra strings are of the type q s — q s , 
with sea quarks and antiquarks at their ends (Fig. |3|). These configurations correspond to 
multiple inelastic scattering in the S'-matrix approach, the number of strings being equal to 
twice the number of inelastic collisions. The contribution of each configuration to the cross 
section is determined using the generalized eikonal approach (see below) or the perturbative 
Reggeon calculus |5C] in hadron-hadron collisions and the Glauber- Gribov model [ |5l"| |52^ ] in 
collisions involving nuclei. 



For A-A collisions, the rapidity distribution of secondaries is given by [28 
dN AA 



dy 



n A 



j^{qq)v v -qv ' ( y \ + N <iv 



(y) 



+ 2 (n -n A ) N qs - q °(y), 



(1) 



6 The value of 8000 charged particles per unit rapidity at y = was the preferred value for many models 
before 1995. Indeed only the SFM gave values close to 2500. Now, by different although probably related 
reasons, several models have lowered their predictions to values close to the SFM one. 
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where N(y) are the rapidity distributions of produced particles in the individual strings 
stretched between the projectile (p) and target (t) nuclei, n A is the average number of 
wounded nucleons of A and n is the average number of nucleon-nucleon collisions. Both 
n\ and n are computed |5^] in the Glauber model. For instance, for minimum bias collisions 

A 2 a NN A 4 / 3 



_ AVnn 
n = 



(2) 



a AA 7r(2i? A ) 2 4 ' 
with ctnn and <jaa the nucleon-nucleon and nucleus-nucleus cross sections respectively; for 
central collisions, 



ctaa ^ tt-RI 



n 



A 4 / 3 . 



(3) 



If all strings would have the same plateau height (i.e. the same value ofiV(O)), dN AA /dy\ y=0 
would increase like A 4//3 . However at present energies the plateau height of the q s — q s strings 
is much smaller than that of the (qq) v — qv ones, and the first term in (|l]) dominates. One 
obtain in this way the result of the Wounded Nucleon Model (WNM) |54|| . At SPS energies 



only for central collisions some departure of the law [29] 

dN AA 



dy 



OC 71 A 



(4) 



Ill- 



is expected, and indeed has been seen in the experimental data [ 

At higher energies the contribution of the sea strings becomes increasingly important, 
not only because their plateau height gets higher but also due to the need to introduce 
multistring configuration in each nucleon-nucleon collision. If the average number of strings 
in each nucleon-nucleon collision is 2k (this number can be computed in the generalized 
eikonal model), the total number of strings is 2kn and (|l|) is changed into 



dN AA 
dy 



N {m)t p -q> ( y ) + fftfr-iOQit* (y) + (2k - 2)N q °- q s (y) 



+ 2 k [n-n A ) N qa - q °(y). 



(5) 



The hadronic spectra of the individual strings N(y) is obtained from a convolution of 



momentum distribution functions and fragmentation functions [27|. Both functions can be 
determined to a large extent from known Regge trajectories. 

For RHIC and LHC energies k ~ 2 and 3 respectively. Using these values in (||) it is 



obtained [57|: 
dN ss 



dy 
dN ss 



dy 



y=0 



y=0 



170, 
500, 



dN 



PbPb 



dy 
dN PbPh 



dy 



y=o 



y=o 



1890 at = 200 GeV per nucleon, 

7900 at \fs = 7 TeV per nucleon, (6) 
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for charged particles in central (tTa > 28 in S-S and 200 in Pb-Pb, corresponding to b ~ 0) 
A- A collisions (in |jM] a value of 8500 for Pb-Pb, b < 3 fm, at \/s = 6 TeV per nucleon is 
given) . 

In these results no semihard collisions were taken into account. The inclusion of this 
kind of collisions cannot modify significantly the numbers in (|6|) since the total number of 
strings is constrained by unitarity. The fact that some of the q s — q s strings can be the 
result of a semihard gluon-gluon interaction, will affect the p± distribution of the produced 
particles. However, average multiplicities are practically unchanged if one neglects changes 
from energy-momentum conservation due to the larger p± in the semihard contribution. 



Other effect not taken into account in these first estimations done in 1991 57] is shadowing 
corrections, which can be of importance at RHIC and LHC energies. The physical origin of 
shadowing corrections]] can be traced back to the difference between the space-time picture 
of the interaction in the Glauber model |5lJ and in Glauber-Gribov field theory [|52| ]. In 
Glauber we have successive collisions of the incident hadron to explain multiple scattering 
in hadron- nucleus interactions, while in Gribov theory simultaneous collisions of different 
projectile constituents with nucleons in the target nuclei are considered. Nevertheless, the h- 
A scattering amplitude can be written as a sum of multiple scattering diagrams with elastic 
intermediate states, which have the same expressions in both cases. In addition to these 
diagrams there are other ones which contain, as intermediate states, all possible diffractive 
excitations of the projectile hadron, whose influence at SPS energies is small. The size of the 
high mass excitations of the initial hadron is controlled by the triple Pomeron coupling. The 
value of this coupling, determined from soft diffraction experimental data, allows to describe 
hard diffraction measured at the Hadron-Electron Ring Accelerator (HERA) at DESY and 
also the size of the shadowing effects in the nuclear structure functions at small x |^j|] . These 
considerations imply a reduction of particle densities at midrapidity of a factor 2 at RHIC 
and 3 at the LHC @. 

This shadowing can be alternatively seen as a way of introducing the interaction among 
strings, see next Section. In (]!]) and (||) it is assumed that strings fragment independently 
one from each other. As the number of strings grows with the energy of the collision, with the 
size of the projectile or the target and with the degree of centrality of the collision, interaction 
of strings is expected at very high energies in central heavy ion collisions. This approach is 
equivalent to take into account the triple Pomeron coupling, whose effects on dN/dy are very 
small at SPS energies and become large at RHIC and LHC. 

In order to include the hard part in the DPM, the eikonal depending on impact parameter 
7 For a discussion on the relation between unitarity, parton saturation and shadowing see for example MM . 
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b and energy is divided in a sum of soft plus hard pieces |6(j, [HJ , 

X (b 2 ,s) = Xs (b 2 ,s)+Xh{b 2 ,s), 

normalized to the corresponding elementary cross sections, 

J d 2 b2 X i{b 2 ,s) = a° , i = s,h; 

in terms of the eikonal, the inelastic cross section for the collision is 



The soft eikonal is parametrized as 

Xs(b 2 ,s) a ° 
and the hard one as 



d z b 1 - e 



-2x(bV) 



8tt [c + a' log (s/so)] 

2 ^ 

Xh(b ,s. 



exp 



4 [c + a' log (s/s Q )} 



8^d eXP {-4d) ■ 

The soft input is a soft Pomeron with a linear trajectory, a s (t) = 1 + A s + a't, 



(7) 



(8) 



(9) 



(10) 



(12) 



and the hard cross section cr? is calculated from PQCD using a lower p± cut-off and con- 
ventional structure functions. Unitarity of the cross section is explicit in (0), which can be 
expanded as 

a in = Jd 2 b cr(l c ,m c ,b 2 ,s), (13) 

i c +m c >l 

the sum running over l c soft elementary collisions and m c hard ones. 

DPMJET is a Monte Carlo code for sampling hadron-hadron, hadr on- nucleus, nucleus- 
nucleus, lepton-hadron and lepton-nucleus collisions at accelerator and cosmic ray energies 
[30]. It uses the DPM for hadronic and nuclear interactions, the hard part being simulated 



using PYTHIA [62] and, in its latest version DPMJET-II.5 [31], one of the most recent sets 



of parton distribution functions, GRV-LO-98 63]. The code includes intranuclear cascade 
processes of the created secondaries with formation time considerations, and also nuclear 
evaporation and fragmentation of the residual nucleus. 

In the first versions of the code, in addition to diagrams where the valence diquarks at 
the end of one string fragment into hadrons preserving the diquark, diquark breaking was 
allowed. This is the so-called popcorn mechanism, see Fig. |]. However the mechanism is 
not enough to explain the large baryon stopping observed in A-B collisions. For this reason 
in the DPM new diagrams [19, ^TJ, 64, |65| for diquark breaking, like the one in Fig. ||, have 
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been proposed and discussed. In both Figs. ^| and [5] the dashed line is the string junction: 
at large N c a baryon can be pictured |3(| as made out of three valence quarks together with 
three strings which join in the string junction^] . These diagrams, included in DPMJET-II.5, 
shift the baryon spectrum to the central rapidity region and also produce an enhancement of 
strange baryons and antibaryons. 

In the code the presence of diquarks and antidiquarks at sea string ends is also included. 
This increases baryon and antibaryon rapidity densities and, due to energy-momentum con- 
servation, reduces that of pions. 

The results of DPMJET-II.5 for charged particles in central Pb-Pb collisions at RHIC (3 
% more central events) and LHC (4 % more central events) are [31] 

RHIC m 



dN 
dy 



= 1280, 

y =o d y 



LHC 

= 2800. (14) 

y=0 



These values agree with the ones computed in the DPM [58], see above. The previous version 
of the code [|(| gives a higher value at LHC, dN c h/ dr)\ v= Q = 3700 |^] (although for yfs = 6 
TeV per nucleon and b < 3 fm). This reduction is due to the inclusion of new diagrams and 
to the energy-momentum conservation consequences of the inclusion of (qq) s — (qq)s strings. 
Notice than the value obtained in the code is much smaller than that obtained in the DPM 
using (H). This fact is essentially due to energy- momentum conservation, which prevents 
some of the (n,n\,k) configurations to take place. 



3 The String Fusion Model 



The SFM [E3, 34, Bq] is based on QGSM, a model which is quite similar to DPM with only 



minor differences. The main ingredient added in the SFM is the fusion of strings [68, 66]. 
The basic idea is that strings fuse as soon as their transverse position come within a certain 
interaction area, of the order of the string proper transverse dimension as dictated by its 
mean p±. In a Monte Carlo approach, such a picture can be realized by assuming that 
strings fuse as soon as partons which act as their sources, have their transverse positions 
close enough. In this language the fusion probability is determined by the parton transverse 
dimension, that is, by the parton-parton cross section. Energy conservation can be taken 
into account by distributing the available energy among these active partons, as it has always 
been done in string models |27], [32|. Then the emerging strings occupy different intervals in 
rapidity space, determined by the energy-momentum of their sources. The fusion of strings 
may only take place when their rapidity intervals overlap. In particular, for two pairs of 



8 Using some supergravity solution and the recently conjectured duality between gauge and string theory 
the large N c baryon wave function has been constructed from N c strings connected via a junction p^ . 
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partons from the projectile and target with rapidities y%, 1/2 and y[, y' 2 respectively, the two 
corresponding strings fuse in the interval [max{y^, y' 2 }, minjyi, 2/2}]- If this interval becomes 
small the resulting object will have will have its total energy of the order of a typical hadron 
mass and, as with ordinary strings, is no more a string but rather an observed hadron. The 
exact value of the minimal string energy and thus of its minimal rapidity length is taken the 
same as for ordinary strings. 

The color and flavor properties of the formed strings follow from the properties of their 
ancestor strings. The fusion of several quark-antiquark q — q strings produces &Q — Q complex 
with color Q (quadratic Casimir operator of the representation Q 2 ), which is determined by 
the SU(3) color composition laws. For example, the fusion of two q — q triplet strings produces 
a [3] string (that is, a diquark-antiquark string) with probability 1/3 and a [6] string with 
probability 2/3 ([3] <g> [3] = [6] © [3]). On the other hand, if two triplet strings with opposite 
color flux directions fuse (a quark [3] state fuses with a [3] antiquark state), either colorless 
states at the end of the new string or a [8] string are formed with probabilities 1/9 and 8/9 
respectively ([3] © [3] = [1] © [8]). The flavor of the fused string ends is evidently composed 
of the flavor of the partons sitting there. As a result of string fusion, we thus obtain strings 
with arbitrarily large color and differently flavored ends, in accordance with the probability 
to create the color Q from several (anti)quarks. Crude characteristics of hadron interactions 
depend only on the fact that the total number of strings of whatever color in a given transverse 
area becomes limited because of their fusion. In other words, string density cannot grow 



infinitely but is bounded from above [38|. More detailed properties of hadron spectra require 
knowledge of a particular manner in which the new fused strings decay into hadrons. 

In all color string models, it is assumed that the homogeneous color field corresponding 
to the strings creates pairs of colored partons, which neutralize this field and provide for 
its subsequent decay. The basic formula which describes the probability of such a process 
is taken in the spirit of the famous Schwinger expression for the probability to create an 
electron-positron pair in a constant electromagnetic field |7l|, (72|. With a constant color 
field which originates from two opposite color charges Q and —Q (8-vectors in SU(3)), the 
probability rate to create a pair of partons with color charges C and — C, flavor / and mass 
Mf for unit string length is assumed to be given by 

~kr K At {kQC) exp r w J ' ( } 

The parameter At has the meaning of the string transverse area and M± is the transverse 



mass, k is proportional to the string tension k, 

irkQ 2 



(16) 
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In |7^, ^] the strings of high color, denoted as color ropes, break as a result of successive 
production of qq pairs, which gradually neutralize the color flux of the string until it breaks. 
In SFM, however, the process considered is that of creation of a pair of parton complexes with 
color Q equal to that of the ends of the string. This is the main contribution to the breaking 
of the string for low values of Q [^] . As, in the Monte Carlo code, fusion of strings is taken 
into account in an effective way and only fusion of two strings is allowed, the mechanism of 
string breaking for high color strings is the one just mentioned. It is also assumed that [3], 
[3] j [6] and [8] strings have the same transverse area[] . The string tension is proportional to 

Qf3 ]= 4/3 = Q|, Qf 8] =3, Qf 6] = 10/3. (17) 

So, approximately K[ 8 j ~ K[ 6 j ~ 2.5 ft[ 3 ] = 2.5 k^. 

As can be inferred from what has been presented above, fusion of strings leads to an 
enhancement of baryon and antibaryon production, due to the possibility of having (qq) and 
(qq) at the end of the fused strings. In DPM a similar mechanism is introduced considering 
the possibility of diquarks in the sea, as mentioned in the previous Section. In addition to this 
mechanism there is another source of baryon enhancement, due to the larger tension of fused 
strings (17) which, through ( |l6| ) and (|l5|), implies a more efficient production of heavy quarks 
and diquarks, and of higher p±. Therefore, it is also expected heavy flavor enhancement and 
some increase of transverse momentum. 

Another important consequences of string fusion is the possibility of producing particles 
in collisions involving nuclei, outside the nucleon-nucleon kinematical region, the so-called 
cumulative effect (part of these effect is usually addressed to the Fermi motion of nucleons 
inside nuclei). In fact the resulting fused string has an energy- momentum corresponding to 
the sum of the energy-momenta of its ancestor strings, which can be larger than the energy- 
momentum available in an isolated nucleon-nucleon collision J77|| . 

In the SFM code, the nuclear parton wave function is taken as a convolution of the parton 
distribution in a nucleon with the nucleon distribution in a nucleus. In this way, hadrons and 
nuclei are treated in a similar way, different from what DPMJET does. Also, in a previous 



version of SFM [23] most of the computations where done at SPS energies and no hard part 
was considered. This part is now introduced in the code in a standard way [35| and modifies 
the central rapidity region at energies higher than those of SPS. 

The probability of fusion of two strings is controlled by the parton-parton cross section, 



2vrr 2 . (18) 



9 This assumption is a very strong one; other possibilities will be discussed in Section 8 in relation to 



percolation of strings QQ^Q. 
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Its numerical value is fixed to a v ~ 8 mb, in order to reproduce the A enhancement seen 
in central S-S and S-Ag collisions at SPS energies J78|. This value, which means r ~ 0.36 
fm, has been obtained implementing in the code fusion of only two strings and therefore has 
to be considered as an effective one. The actual transverse size of a string should be less, a 
more realistic one being r ~ 0.2-^0.25 fm |l(]], a value which agrees with other considerations 
[71 [f|,|(J. 

For the purpose of this review, the most important consequence of fusion of strings is 
that it suppresses total multiplicities, reducing the number of pions in the central rapidity 
region, although the rapidity distribution becomes larger at the extreme of the fragmentation 
regions. In the predictions done with the previous version of the model p3| , no hard part 
was included and the charged particle densities at midrapidity for central (b = 0) Au-Au 
collisions were 



dN 
dy 
dN 
dy 



y=0 



1000 at y/a = 200 GeV per nucleon, 

1900 at = 6.3 TeV per nucleon. (19) 



y=o 

The corresponding values without considering fusion of strings were 1850 and 4000 respec- 
tively. A strong suppression of the central density is produced (note the agreement with 
the values from DPMJET). Including the hard part ||3Ef| , the values for collisions of charged 
particle densities corresponding to the 5% more central events^ are 

dN AuAu 

= 910 at y/s = 200 GeV per nucleon, 



dy 
dN PhPh 



dy 



= 3140 at = 5.5 TeV per nucleon, (20) 

y=o 



and now the suppression due to string fusion is smaller (the corresponding values without 
fusion are 1300 and 3690 respectively). The reason for this is that the strings coming from 
hard scatterings do not fuse in the code. At LHC a large proportion of strings are hard ones 
and therefore the relative size of suppression is smaller. The hard strings have a size ~ l/p± 
and indeed should interact and fuse, although with smaller probability than the soft ones. 
Effects of overlapping of strings will be further discussed in Section 8. 



4 The Relativistic Quantum Molecular Dynamics Model 



RQMD H, |7|, H| is a semiclassical microscopic approach which combines classical prop- 
agation with stochastic interactions. Strings and resonances can be excited in elementary 
10 In the model, this translates into b < 3.2 fm for Au-Au at RHIC and 6 < 3.3 fm for Pb-Pb at LHC. 
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collisions, their fragmentation and decay leading to the production of particles. The nature 
of the active degrees of freedom in RQMD depends on the relevant length and time scales of 
the processes considered. In low energy collisions (around 1 GeV per nucleon in the center of 
mass) RQMD reduces to solving transport equations for a system of nucleons, other hadrons 
and eventually resonances interacting in binary collisions and via mean fields. At large beam 
energies (> 10 GeV per nucleon in the center of mass) the description of a projectile hadron 
interacting in a medium (a cold nucleus) as a sequence of separated hadron or resonance col- 
lisions breaks down. A multiple collision series is formulated on the partonic level, following 
the Glauber-Gribov picture. In RQMD these multiple collisions correspond to strings formed 
between partons of the projectile and target, including sea quarks and antiquarks. The string 
excitation law dP oc dx + /x + is the same used in FRITIOF [33]. The decay of elementary 
color strings is done using JETSET [|81[J . Rescattering is included: four classes of binary 
interactions, BB, BM, MM and BB (B denoting baryon, M denoting meson) are considered. 

One of the main ingredients of RQMD is the inclusion of interaction of strings by means 
of formation of color ropes, see previous Section, when there are overlapping strings. These 
ropes are chromoelectric flux tubes whose sources are charge states in representations of color 
SU(3) with dimension higher than the triplet one. They are equivalent to the fused strings of 
SFM. As already mention, as a simplification in SFM only fusion of two strings is considered, 
as an effective way to take into account string interaction. In RQMD all possibilities are 
considered. The breaking of these higher color strings proceeds through successive production 



of qq pairs [72, 73] due to the Schwinger mechanism. 

As in the case of SFM, introduction of color ropes in RQMD leads to heavy flavor and 
baryon and antibaryon enhancement. In version RQMD 2.3 the model reproduces the SPS 
rapidity distributions of h _ , K°, A, A, H~ and H + . It slightly underestimates the yields of 
£l~ and + (by less than a factor 2). Also it is able to reproduce the m± spectrum of all 
these particles. Let us mention that independent string models are not able to reproduce 
these slopes. 

The formation of color ropes leads to a strong suppression of central rapidity distributions. 
The prediction of RQMD for central (6 = 3 fm) Pb-Pb collisions at RHIC is dN/dy\ y=0 ~ 700 



[26, |82[| . This number is lower than the value of SFM, 910. The reason for that probably has 
to do with the strong fusion probability used in RQMD. The effect of this strong string inter- 
action in some observables (like antibaryon enhancement) is compensated by other processes 
(a large BB annihilation). 
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5 The Heavy-Ion Jet Interaction Generator (HIJING) 

In HIJING [0] the soft contribution is modeled by diquark-quark strings with gluon kinks 
induced by soft gluon radiation, in a way very similar to the FRITIOF model @|. Since this 
model treats explicitly minijet physics through PQCD, the transverse momentum in string 
kinks due to soft processes is limited from above by a minijet scale p±o = 2 GeV/c. Gluon 
radiation is extended to the hard part of high p±, which, together with the use of momentum 
distribution functions for partons similar to those of DPM, constitutes a difference with 
FRITIOF. Strings decay independently by means of the JETSET pl| ] routines. In addition 
to the low p±_ < pj_q gluon kinks, HIJING includes an extra low p± transfer between the 
constituent quarks and diquarks at the string ends. This extra p± is chosen to ensure an 
smooth extrapolation in the p± distributions from the soft to the hard regime. 

Multiple minijet production with initial and final state radiation is included along the 
lines of the PYTHIA model |32|. First, the cross section for hard parton scattering aj e t is 
computed in PQCD at leading order (LO), using a K- factor ~ 2 to simulate higher order 
corrections. The eikonal formalism, see Section 2, is employed to calculate the number of 
minijets per inelastic nucleon-nucleon collision. For A-A collisions at impact parameter b the 
total number of jets is given by 



atAA /ia A T AA (b) 



with 



T AA (b) = J d 2 b' T A {b - b')T A (b'), (22) 
T A (b) being the nuclear profile function normalized to 1 and a AA (b) ~ J d 2 b {1— exp [— a^A 2 T A 



the A-A cross section |29| , 53] for impact parameter b. For central collisions 6 = 0, cr AA (b = 
0) ~ 1 and 

N^(b) « °Jet oc A 4//3 . (23) 

Therefore, at high energies and for central nucleus-nucleus collisions there will be many 
minijets. 

In the model, jet quenching Q is included to enable the study of the dependence of 
moderate and high p± observables, on an assumed energy loss per unit length dE/dx of high 
energy partons traversing the dense matter produced in the collision. The effect of including 
jet quenching is a moderate enhancement of particle production in the central rapidity region 
and to diminish the yield in the fragmentation regions. Furthermore, in the last version of 
the model the mechanism of string junction migration explained in Section 2 is included, 
in order to shift baryons from fragmentation to central rapidity regions. 
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The results for charged densities at midrapidity in central (b < 3 fm) Au-Au collisions at 
= 200 GeV per nucleon are shown in Fig. ||. The different curves refer to different versions 
of the model with and without quenching and shadowing of the nucleon structure functions in 
the nucleus pi M. For LHC, HIJING predictions @ lie in the range 5000-^7500 depending 



on structure functions used and quenching included or not. 

6 Perturbative Quantum Chromodynamics and Hydrodynam- 
ical models 

It has been argued [46] that the initial state (the initial distribution of partons) in a high 



energy heavy ion collision could be computed using PQCD. Several groups ]44] , f45| have 
developed models along this line. Concretely, Eskola et al. have computed [H] charged 
densities and transverse energies al midrapidities, using PQCD at some given scale which is 
taken to be equal to a saturation scale, the scale at which parton distributions stop their 
increasing at small x. 

In ultrarelativistic heavy ion collisions the number of produced gluons and quarks with p± 
greater than some cut-off po, Naa (b,Po, y/s), increases when po decreases, when the size of 
the nuclei increases or b decreases, see (^) and (23), and when ^/s increases due to the small 



x enhancement of parton distribution functions. Shadowing of nucleon structure functions in 
nuclei will decrease ./Vaa (b,po, y/s), but next-to- leading order (NLO) corrections will increase 
it. At sufficient large cut-off po 3> Aqcd the system of produced gluons is dilute and usual 
perturbation theory is applicable. However, at some transverse momentum po = p sa ± the 
gluon and quark phase space density saturate and no further increase is expected. In this 
case one may conjecture that evaluation of the number of charged particles N c h and transverse 
energy Et using QCD formulae at this saturation scale p sa t gives a good estimate of the total 
N c h and Et (partons with p± 3> p sa t are rare, partons with pj_ <C p sa t saturate and contribute 
little to the total Et)- 

In [44], first Naa (b = 0,po, y^) for |y| < 0.5 is computed using standard PQCD ex- 
pressions at LO. Nuclear effects on parton distribution functions are implemented using the 



EKS98 parameterization [84] of nuclear corrections. To simulate NLO contributions, a K- 
factor K = 2 is used. The scale in the PQCD calculation is fixed from considering that at 
saturation Naa(^ = 0,p sa t, partons, each one with transverse area ^/p^ a t, fill the whole 
transverse area ttR\, 

Naa (6 = 0, p sa t , <Js) = p 2 sat R\ . (24) 

In Fig. | A AA (b = 0, p , y/s) is plotted for A = 208 as a function of po at SPS, RHIC and LHC 
energies. The dashed curve is PqR\- The intersection points give us p sa t at the corresponding 
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energies. Of course, all this is only valid as long as po 3> Aqcd for perturbation theory to 
be justifiedly used, which is doubtful at SPS and RHIC (see the Figure; the saturation 
momentum are ~ 0.5, ~ 1.4 and ~ 2.3 GeV/c at SPS, RHIC and LHC energies respectively). 
The values of Ni = Naa (b = 0,p sa t, \/s) and p sat can be well fitted by the expressions 

Ni = 1.383 A ' 922 (a/s) ' 383 , (25) 
p sat = 0.208 A ' 128 (^i) - 191 GeV/c. (26) 

The initial state computed in this way very nearly fulfills the kinetic thermalization condition 
for bosons, e/n = 2.7 T (the number of gluons is much larger than that of quarks), and there 
also is some justification to consider that further hydrodynamical expansion is locally thermal, 
i.e. entropy conserving. Thus initially the entropy Si = 3.6 Ni (ideal system of bosons). For 
the final hadronic gas Si = Sf — 4 Nf so that Nf = 0.9 Ni, i.e. the number of hadrons in 
the final state is, up to 10 % corrections, equal to the number of initially produced gluons at 



the scale p sa t- The multiplicity prediction []44j, 26] 



N ch = | 0.9 Ni (27) 

is directly obtained from (^) and plotted as a dashed line in Fig. ||. The values at RHIC 
and LHC for central (b = 0) Pb-Pb collisions are 900 and 3100, not very different from those 
obtained by DPM, DPMJET and SFM on very different grounds^ . 

7 Other models 

In this Section we would like to comments on some other models. The fact that the rest of 
the models are joined together, do not mean at all that they are less important or successful 
than the mentioned ones. It is simply the shortage of space which prevents us from a longer 
study. 



The VENUS model [24] is an extension of DPM. The main difference is the inclusion of 
diagrams in which there is two color exchanges, the first one providing two (qq) v — q v strings, 
one of the being intermediate, because a second color exchange breaks the diquark, giving a 
different {qq) v — qv string and a double string which consists in a forward moving quark linked 
to two backward moving quarks, see Fig. ^. These diagrams become increasingly important 
with a growing number of inelastic collisions, as in h-A or A-B, although of course they are 
also present in N-N collisions. They enhance stopping power, shifting the baryon spectrum 
towards midrapidities. 



11 Theoretical models based on a semiclassical treatment of gluon radiation ]8^] by partons in the colliding 
nuclei give values compatible with these ones, see pq]. 
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VENUS gives large values for central rapidity densities. At = 6 TeV per nucleon 



for central (6 < 3 fm) Pb-Pb collisions its result is dN c h/ '^77^=0 = 8400 [22]. The model 
has lately been extended to deal with 7*-7*, 7*-h, u-h, h-h, h-A and A-B collisions in the 
same unified approach j|l] . An unique Pomeron describes both soft and hard interactions by 
means of the evolution of structure functions from some properly chosen initial conditions. 



The new model [41], denoted by NEXUS, has not given values for LHC yet. Preliminary 
predictions for central (b < 2 fm) Au-Au collisions at RHIC are dN c h/dy\ y= o ~ 1100 [E(J. In 
this model particles and resonances produced in string fragmentation are allowed to rescatter 
and, if more than two of them are close enough, joined into a quark cluster which is decayed 
isotropically p7|]. 

The LUCIAE event generator E3] (Lund University and China Institute of Atomic En- 
ergy) is a version of the FRITIOF model |33| where collective interactions among strings 
and rescattering of the produced particles are included. The collective interactions are in- 
corporated following |pg| l. LUCIAE has studied several observables at SPS comparing to 
experimental data, but, to our knowledge, has not worked out predictions for rapidity densi- 
ties at RHIC and LHC. In any case, the inclusion of collective string effects produces very fast 
particles (7j| at the extreme of the phase space, similar to the cumulative effect which occurs 
in the SFM. Also, simply by energy-momentum conservation, a suppression of particles in 
the central rapidity region should happen. 

The Ultrarelativistic Quantum Molecular Dynamics (UrQMD) [39, ^] is a microscopic 
hadronic approach based on the covariant propagation of mesonic and baryonic degrees of 
freedom. It allows for formation of strings and resonances, and rescattering among them and 
of the produced particles. In this aspect it is quite similar to RQMD. In the low energy region, 
i.e. i/s < 2 GeV per nucleon, the inelastic cross sections are dominated by s-channel formation 
of resonances, which decay into particles isotropically in their local rest frame according to 
their lifetimes. A large variety of baryonic and mesonic states have been incorporated in 
the model. All corresponding antiparticles are included and treated on the same footing. At 
higher energies strings are considered. Much attention is paid in the model to the intermediate 
energy region between AGS, the Alternating Gradient Synchrotron at BNL (y/s ~ 5 GeV 
per nucleon), and SPS, in order to achieve a smooth transition between the low and high 
energy regimes. The prediction of the model [EH, [2(| for central (b < 3 fm) Au-Au collisions 



at RHIC is shown in Fig. 1C. The central density of charged pions is ~ 750, being for all 
charged particles ~ 1100. 

Another approach using also UrQMD is the model denoted by VNI+UrQMD, where a 
combined microscopic partonic/hadronic transport scenario is introduced p9| . The initial 
high density partonic phase of the heavy ion reaction is calculated in the framework of the 
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parton cascade model VNI f90|| , using cross section obtained from PQCD at LO (see |9lJ] 
for a discussion of the uncertainties introduced by NLO corrections and K-factors in VNI 
and other parton cascade models). The partonic state is then hadronized via a configura- 
tion space coalescence and cluster hadronization model, and used as initial condition for a 
hadronic transport calculation using UrQMD. In Fig. Ill] the time evolution of parton and 



on-shell hadron rapidity densities for central (b < 1 fm) at RHIC can be seen [89, 26]. From 
this curve, the charged particle rapidity density is ~ 1000. When, instead of calculating the 
initial phase using a parton cascade approach, QGP formation is assumed, this plasma is 
evolved hydrodynamicalfy until hadronization, and then UrQMD is used for hadronic trans- 
port (the model considers a first order phase transition and is denoted as Hydro+UrQMD 
p2fl ), a smaller charged particle central rapidity density is obtained for central (b = 0) Au-Au 
collisions at RHIC, ~ 750. 

Using also VNI there is the model denoted by VNI+HSD |93|], where HSD stands for 



Hadron String Dynamics model [94]. This model involves quarks, diquarks, antiquarks, 
antidiquarks, strings and hadrons as degrees of freedom. The parton cascade model VNI is 
extended by the hadronic rescattering as described by HSD. Its results for central (b < 2 
fm) Au-Au collisions at RHIC are shown in Fig. where VNI and HSD predictions are 
the four plots at the top and VNI+HSD are the two plots at the bottom. It is worth noting 
that VNI+HSD gives almost the same results as HSD, which already includes final state 
interactions. This fact, also observed in DPM and SFM (see Section 2) is based on unitarity, 
which controls the number of inelastic collisions independently of their soft or hard origin. 
The total number of charged particles at midrapidity is ~ 1150. 

Another RHIC prediction []95| , p6f| comes from a modification of the HIJING model to 
include a parton cascade model [Bq] and final state interactions based on the ART model 



[97]. In this model (HIJING+ZPC+ART) the central rapidity density of charged particles 
for central (b = 0) collisions at RHIC is of the order 1100. 

The assumption that local thermodynamical equilibrium is attained by the system of 
two heavy ions colliding at high energies is a basic hypothesis of macroscopic statistical and 



thermodynamical models |9S, 99, 10C, |101[| (see |102[| for a discussion on statistical equilibrium 



in heavy ion collisions). This idea comes from a long time ago [103, 104, |105| [. Following 
||, [L01| , |26|| the statistical model treats the system as a grand canonical ensemble with two 
free parameters, a temperature T and a chemical potential Interactions of the produced 
particles are taken into account considering a excluded volume correction, corresponding 
to repulsion setting in for all hadrons at a radius of 0.3 fm (a hard core). Hadron yield 
ratios resulting from this model are in reasonable agreement with SPS central Pb-Pb data. 
Taking these results and looking at the expected phase boundary between the QGP and the 
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hadron gas, the hadrochemical freeze-out points are where one expects, this being suggestive 
that hadron yields are frozen at the point when hadronization of the QGP is complete. 
This gives for RHIC a freeze-out temperature of 170 MeV, the same as found at SPS. The 
chemical potential is expected to be small, 10 MeV used as an upper limit. Strangeness 
and I3 conservation then require values of fj, s = 2.5 MeV and fij 3 = —0.2 MeV. In order to 
predict absolute yields one has to estimate the volume per unit rapidity at the time when 
hadronization is complete. Starting from an initial temperature of Tj = 500 MeV at a time 
r = 0.2 fm/c and using a transverse expansion with (5 = 0.16, this volume is 3600 fm 3 at the 
freeze-out temperature of 170 MeV. The number of charged pions per unit rapidity at y = 
is 1260, that of charged kaons 194, that of protons 62 and that of antiprotons 56. Modifying 
the freeze-out temperature to 160 MeV results in a reduction of hadron yields about 10 %. 
For central Pb-Pb collisions at the LHC, taking again Tf = 170 MeV and ^b = 10 MeV 
and performing similar calculations, the fireball volume per unit rapidity at chemical freeze- 
out is 14400 fm 3 , resulting in 5000 charged pions, 770 charged kaons, 250 protons and 220 
antiprotons per unit rapidity. The total charged particle density is 



dN 
dy 



= 7560. (28) 

y=0 



As at RHIC, small modifications of Tf and ^b lead to small changes in this prediction. 
Both at RHIC and LHC the centrality of these results should correspond roughly to the 
centrality of the SPS experimental data which were fitted to extract the parameters used in 
the predictions. As different experiments consider different centrality criteria, this is not fully 
determined, so let us take 5 4- 10 % as an estimate. 



The quark coalescence model [26] assumes that at RHIC a QGP will be produced in 
the collision, which will expand and cool, hadronization proceeding via quark coalescence 



as described by the ALCOR model [106]. In this nonlinear coalescence model, subprocesses 
are not independent, competing one with each other. The coalescence equations relate the 
number of hadrons of a given type to the product of the numbers of different quarks from 
which the hadron consists. The main predictions of the model relates different particle ratios, 
but unfortunately the absolute values cannot be obtained in the model. 

Finally, let us mention an extrapolation done by the WA98 Collaboration |5(| from central 
pseudorapidity densities measured at SPS at different centralities. The predicted maximum 
charged pseudorapidity density for central Pb-Pb collisions is 1000 at RHIC and 2500 at 
LHC, for ~ 10 % more central events. 
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8 Percolation 

In many models, multiparticle production at high energies is described in terms of color 
strings stretched between the projectile and target, see previous Sections. In principle, these 
strings fragment independently, the only correlation among them being energy-momentum 
conservation. However, with growing energy, centrality and/or size of the colliding particles, 
the number of strings grows and one expects that the hypothesis of independent fragmentation 
is no longer valid, interaction among them becoming essential. For these reasons we have 
seen in Section 3 and 4 different ways of taking into account such interaction. In particular, 
in SFM or RQMD the strings fuse or form color ropes, in such a way that the transverse size 
of the new string or rope is the same as that of the original strings. In this case it can be 
shown that there is no phase transition. However, other possibilities could be discussed. It 
could be the case that the new string would have a transverse size corresponding to the sum 
of the sizes of the original strings. In this case, a first order phase transition occurs |75[| . 

An alternative and natural way to the formation of a (non-thermal) QGP is percolation 
of strings [16, 107], which in any case can be used as an estimation for the failure of the 



independent fragmentation hypothesis. This a purely classical, geometrical mechanism. At 
a given energy and impact parameter b in an A-B collision, there is an available transverse 
area. For simplicity, let us take A=B and b = 0, so this area is ttR\. Inside it, the strings 
formed in the collision can be viewed as circles of radius ro- Some of the circles may overlap 
and then form clusters of strings. Above a critical density of strings, percolation occurs, so 
that clusters of overlapping strings are formed through the whole collision area. Percolation 
gives rise to the formation of a collective state, which can be identified as QGP, at a nuclear 
scale. The phenomenon of continuum percolation is well known and has been used to explain 
many different physical processes |108fl . The percolation threshold r\ c is related to the critical 
density of strings n c by 

rj c = ttvq n c , (29) 
n c = —4j at 6 = 0, (30) 

7T it a 



A 



sec 



where iV c is the number of exchanged strings and ro — 0.2 -j- 0.25 fm [16, [74|, 79, 
Section 3. r/ c has been computed using Monte Carlo simulation and direct connectedness 



expansions [105]. The results lie in the range 1.12 -j- 1.18 using step functions for the profile 
function of the nucleus (i.e. strings homogeneously distributed in the whole transverse area 
available). The use of Woods-Saxon or Gaussian nuclear densities leads to higher values of 
rj c [ |109| ], up to ~ 1.5. The corresponding value of n c lies in the range 6 -f- 12 strings/fm 2 . 
In Table [l] we show the number of strings exchanged in different central (6 = 0) collisions 
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together with their densities, as obtained in the SFM |]35|| . It is seen that percolation could 
already occur for Pb-Pb at SPS. At RHIC and LHC, even collisions between much lighter 
nuclei could lead to the phase transition. 

Table 1: Number of strings (upper numbers) and their densities (in strings/fm 2 , lower num- 
bers) for different central (b = 0) collisions at SPS, RHIC and LHC energies, as obtained in 
the SFM g|. 



yfs (GeV per nucleoli) 


P-P 


S-S 


Pb-Pb 


19.4 


3.4 


144 


1365 




1.7 


3.4 


9.2 


200 


4.3 


223 


2029 




2.1 


5.2 


13.7 


5500 


5.8 


416 


3469 




2.9 


9.8 


23.4 



Notice that string percolation occurs in two dimensions. Percolation of hadrons was 



proposed long ago [110] as a possible way to reach QGP. However, in this case percolation is 



three-dimensional and the critical density is below even normal nuclear matter density. This 



is in agreement with lattice studies [80] which show that the so-called energy radius of the 
hadron is about 0.2 fm. Therefore the color fields inside hadrons occupy only a few percent 
of the transverse area, (^o/iih) 2 ra (1/5) 2 - This also explains the relative weak string-string 
interaction (for instance the triple Pomeron coupling of Glauber-Gribov theory used in DPM 
and DPMJET, see Section 2). 



Percolation is a second order phase transition [105]. The corresponding scaling law gives 



the behavior of the number of clusters of n strings, u n , in terms of rj, 

(u n ) = n- T T {rf [77 - r) c \) , \n - rfcl « 1, n » 1, (31) 

where r = 187/91 and a = 36/91. The fraction eft of the total surface occupied by strings is 
determined by 

= 1- e - ". (32) 



It can be seen |75j] that the multiplicity /x n due to a cluster of n overlapping strings, compared 
to the multiplicity of one string fii, is given by 

fin = n [i x F(rj), (33) 
1 1 - e -v 

F(n) = J^ r - (34) 
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Prom Table |l] and taking tq = 0.25 fm, the values of F(rj) for central (b = 0) Pb-Pb collisions 
at RHIC and LHC are 0.59 and 0.46 respectively, quite close to the values of the reduction 
factor of multiplicities in DPM due to triple Pomeron couplings, 1/2 and 1/3 respectively, 
see Section 2. A naive calculation for central rapidity densities of charged particles can be 



done multiplying the corresponding values obtained in the SFM [35] without fusion of strings 
for Pb-Pb at b = 0, by these reduction factors F(jf). The results for RHIC and LHC are 910 
and 1980 respectively (380 at y^s = 19.4 GeV with F{jf) = 0.68) and clearly, from the way 
they were obtained, should be considered as lower bounds. 

Percolation, in addition to reduce central rapidity multiplicities and enhance heavy flavor 
production (as string fusion does), also gives rise to other consequences in long range cor- 
relations and transverse momentum correlations [|r(J, and 3/tp suppression [74, 10S]. Also, 
as the energy-momentum of the clusters is the sum of the energy-momentum of the original 
strings, a huge cumulative effect is expected. 

The critical point of percolation is the fixed point of a scale transformation (renormal- 
ization group equation) which eliminates short range correlations, surviving only long range 
ones. Close the that point observables should depend only on 77, and not on energy or nuclear 
size separately. 

By passing, let us notice that a similar intrinsic scale has been proposed |$f| in the 
small x physics domain, related with saturation of structure functions or minijets, which was 
discussed in Section 6. This quantity is defined by 

a2 _ N AA (p sat ) 

with p sa t the transverse momentum at which saturation starts. When the number of partons 
NAA(Psat), each one with a transverse size of the order 7r/p 2 ai , verifies A 2 7r ~ p 2 sat , partons 
cover the whole nuclear area. Physics should depend only on the value of A 2 . Furthermore, the 
effective action |85|| which describes small x physics should become critical at some fixed point 
of some renormalization group, the correlation functions depending only on critical exponents 
determined by symmetry considerations and dimensionality. Comparing with percolation, 
indeed ^/(tttq) = N/{irR\) is formally A 2 , with the exchange of N soft strings by N A \{p sat ) 
partons. Let us indicate that from the arguments exposed above, overlapping partons will 
only cover the whole transverse area ttR\ asymptotically. According to (|32|), the transverse 
area S covered by N A A{Psat) partons of transverse size vr/p 2 at is 

S = ttR\ [l - exp (-A 2 vr/rg)] . (36) 



22 



9 Cosmic Ray Physics and heavy ion accelerators 

Usually it is considered that the highest cosmic ray energies, say 10 15 4 10 20 eV, are much 
higher than energies reached at accelerators. With the advent of RHIC and LHC, this is 
not true any longer. Pb-Pb collisions at RHIC and LHC will reach total energies of ~ 10 15 
and ~ 10 18 eV respectively. This means that, although no participant nuclei larger than Fe 
is expected, there will be collective physics to explore in Cosmic Ray Physics. In particular 
changes in the multiplicity originate changes in the development of atmospheric showers. 
Unfortunately, atmospheric showers are dominated by forward particles, and it is in the 
fragmentation regions where models which in the central rapidity region are quite different, 
are more similar. Nevertheless, collective effects like color rope formation will influence the 
fragmentation regions (for example, the enhancement of the cumulative effect) and have 



observable effects in the development of the shower [77, 112]. Besides, PQCD effects may be 
of importance for the transverse broadening of the shower. It would be convenient to apply 
the different models for multiparticle production to simulations of cosmic ray atmospheric 



showers |113| , |114|1 



10 Conclusions 

In Table || predictions of the different models for charged particle densities produced in central 
Au-Au or Pb-Pb collisions at RHIC and LHC are presented, together with the corresponding 
centrality criteria. Some of the predictions of the models are not available and its place has 
been left empty. In order to estimate the discrepancy between different predictions induced 
by different definitions of centrality, by the fact that some of the results are dN/dy\ y= o and 
other dN / dr/l^o , and by using Au-Au or Pb-Pb at RHIC or slightly different energies at LHC, 
in Table ^ we present results obtained in the SFM [35] for different reactions. From them it 



can be concluded that results should be compared allowing for a 20 4 30 % discrepancy. 



For completeness, we have included predictions from percolation and from the WNM [54]. 
As discussed in Section 8, the former should be considered as lower bounds. The latter have 
been obtained computing in SFM |35|] dN c h/dy\ y= Q for p-p, p-n and n-n collisions at the 
corresponding energies, making an isospin weighted average and multiplying by the number 
of wounded nucleons in a central collision, taking 77a = 200 in Pb-Pb and 77a = 190 in Au-Au 
(which correspond to b ~ 0); the result that we get in this way for Pb-Pb at \/s = 17.3 
GeV per nucleon (SPS) is 370. In the following we will not discuss these two quite naive 
predictions. 

At RHIC energies all the predictions are in the range 700 4 1550 and most of them 
in 1000 4 1100. The lowest value corresponds to RQMD. The reason for that, as already 
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Table 2: Predictions of different model for charged particle densities in central collisions. 
Unless otherwise stated, predictions refer to dN/dy at y = 0, for Au-Au at = 200 GeV 
per nucleon at RHIC and Pb-Pb at yfs = 5.5 TeV per nucleon at LHC. The centrality criteria 
and the references from which the results have been taken are indicated in each case. 



Model 


RHIC 


LHC 


DPM 


1000 


2500 


(n A > 200) © 


(Pb-Pb) 


(*/s = 7 TeV per nucleon) 


DPMJET 


1280 


2800 


© 


(Pb-Pb, 3 %) 


(4%) 


SFM 


910 


3140 


(5 %) © 






RQMD 


700 




©n 


(6 = 3 fm) 




HIJING 


600-1150 


5000 ~ 7500 


(dN/dy at y = 0, b < 3 fm) © §§ 






Eskola e£ 


900 


3100 


(6 = 0) H 


(Pb-Pb) 




HIJING+ZPC+ART 


1100 




©@ 


(6 = 0) 




UrQMD 


1100 




©11 


(6 < 3 fm) 




VNI+UrQMD 


1000 






(6 < 1 fm) 




Hydro+UrQMD 


750 




m 


(6 = 0) 




VNI+HSD 


1150 




© 


(6 < 2 fm) 




VENUS 4.12 




8400 (dN/dy at y = 0, 


© 




6 < 3 fm, i/s = 6 TeV per nucleon) 


NEXUS 


1100 




© 


(6 < 2 fm) 




Statistical 


1550 


7560 


(5-10%) m§ 






WA98 extrapolation 


1000 


2500 


(maximum diV/d?7, ~ 10 %) H 


(Pb-Pb) 




WNM 


560 


1220 




(n A = 190) 


(n A = 200) 


Percolation 


910 


1980 


(6 = 0) 


(Pb-Pb) 
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Table 3: Results obtained in the SFM |}f| for charged particle production in different reac- 
tions, compared with the results for dN/dy at y = in Au-Au at y/s = 200 GeV per nucleon 
(910) and in Pb-Pb at yfs = 5.5 TeV per nucleon (3140), with a degree of centrality of 5 %. 



Reaction and observable 


Result 


Difference (%) 


dN/dy at y = 0, Au-Au at ^/s = 200 GeV 


910 





per nucleon, 5 % of centrality (b < 3.2 frn) 






cW/dr? at = 0, Au-Au at V« = 200 GeV 


730 


-20 


per nucleon, 5 % of centrality (6 < 3.2 fm) 






Maximum dN/dr), Au-Au at ^/s = 200 GeV 


770 


-15 


per nucleon, 5 % of centrality (b < 3.2 fm) 






dN/dy at y = 0, Au-Au at y's = 200 GeV 


1030 


+13 


per nucleon, 6 = 






dN/dy at y = 0, Pb-Pb at = 200 GeV 


960 


+5 


per nucleon, b < 3.2 fm 






dTV/dy at y = 0, Pb-Pb at ^ = 5.5 TeV 


3140 





per nucleon, 5 % of centrality (b < 3.3 fm) 






dN/dij at 77 = 0, Pb-Pb at a/s = 5.5 TeV 


2620 


-17 


per nucleon, 5 % of centrality (b < 3.3 fm) 






Maximum dN/drj, Pb-Pb at ^fs = 5.5 TeV 


2820 


-10 


per nucleon, 5 % of centrality (b < 3.3 fm) 






dN/dy at y = 0, Pb-Pb at V« = 5.5 TeV 


3590 


+14 


per nucleon, 6 = 






dN/dy at y = Q, Pb-Pb at ^/i = 7 TeV 


3330 


+6 


per nucleon, 5 % of centrality (6 < 3.3 fm) 
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commented at the end of Section 4, is the formation of color ropes with a large probability, 
which is required in the model to describe antibaryon enhancement at SPS because a rather 
large baryon-antibaryon annihilation cross section is used. In SFM, which considers a similar 
mechanism (string fusion), the charged density obtained is larger. The low value also obtained 
in Hydro+UrQMD is due to the initial QGP state and the first order transition assumed. 

At LHC the differences among the predictions are larger, more than a factor 3. They lie 
in the range 2500 -f- 8500. Comparing Table ^ with Fig. [I], which summarizes the situation 
before 1996, it can be seen that nowadays predictions tend to gather around the lowest values; 
at the time of p^] only SFM gave predictions below 4000. 

Essentially, models based on parton shower evolution predict larger values. Also statistical 
models obtain very large charged densities. However, the model described in Section 6, which 
uses mainly partonic degrees of freedom, obtains a low value (3100). This is due to the 
saturation of minijet production, which plays the role of an upper cut-off in the number of 
minijets. This saturation is a consequence of unitarity |4q| . 

Unitarity is a basic ingredient that controls the number of soft and hard elementary scat- 
terings (which are no more independent one from each other) in models like DPM, DPMJET, 
SFM and RQMD. In addition to that, energy-momentum conservation reduces the possible 
number of scatterings. Finally, interaction among strings is another collective effect which 
reduces central pion densities. In DPM these interactions are taken into account by means 
of the triple Pomeron. Its coupling is fixed to describe soft diffraction and HERA data. In 
other models as RQMD or SFM the interactions among strings are taken into account via 
the formation of color ropes or fusion of strings. Its strength is essentially fixed to repro- 
duce heavy flavor and antibaryon enhancement at SPS. It is not unexpected that these three 
different forms of quantifying the shadowing give rise to similar predictions for global and 
simple observables as central rapidity densities of charged particles. Probably the knowledge 
of shadowing from small x physics can help to reconcile models based on partonic degrees of 
freedom with those based on strings as degrees of freedom. On the other hand, statistical 
thermal models predict a larger value close to 8000 as a consequence of the rather large vol- 
ume at freeze-out. A reduction in a factor 3 would mean a strong reduction in this volume 
and a large change of the ratio between different particles, or in the temperature and chemical 
potential values. 

In the summary talks of three mayor conferences [||, [|, |115| were quoted 8000, 8000, and 
between 3000 and 8000, for the probable number of charged particles per unit rapidity at 
the center of mass in central Pb-Pb collisions at the LHC. Models with interaction among 
strings obtain a lower value. The interplay between soft and hard physics is one of the main 
issues in the study of strong interactions and, together with the search and characterization 
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of QGP, one of the main goals of Heavy Ion Physics. Doubtless, the new experiments at 
RHIC and LHC will shed light on this subject, even measuring such a simple observable as 
central rapidity densities of charged particles. 
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FIGURES 




Figure 1: Predictions from different models for the charged pseudorapidity density in central 
(b < 3 fm) Pb-Pb collisions at the LHC, taken from (||. 
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Figure 2: Lowest order contribution to p-p scattering at high energies in the DPM. Strings 
are stretched between valence constituents of the protons and hadronize by means of qq pair 
production. 




Figure 3: First higher order contribution to p-p scattering at high energies in the DPM. 
Besides the contribution shown in Fig. ^, now strings stretched between sea constituents of 
the protons appear. 
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Figure 5: One example of proposed diagram for diquark breaking in nucleon- nucleus collisions. 




Figure 6: Predictions for the charged pseudorapidity density in central (b < 3 fm) Au-Au 
collisions at RHIC in the HIJING model, taken from |pq |. 
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Po [GeV] 



Figure 7: 7Vaa(& = 0,po> with pj> > po and |y| < 0.5 for A = 208 as a function of po at 
s/s = 5500 (LHC), 200 (RHIC) and 17 (SPS) GeV per nucleon, taken from Q. The dashed 
curve is PqR\- 
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Figure 8: Rapidity density of charged particles at y = for A = 12, 32, 64, 136 and 208 as 
a function of yfs for central A- A collisions, taken from [|44| . The dotted line is the result of 
(pTp, the solid line is the result of a more rigorous calculation. 




Figure 9: Example of diagram with two color exchanges in VENUS. 
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Figure 10: Predictions for the charged rapidity density in central (b < 3 fm) Au-Au collisions 
at RHIC in the UrQMD model, taken from |§. 
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gure 11: Time evolution of parton and on-shell hadron rapidity densities for central (b < 1 
at RHIC in the VNI+UrQMD model, taken from M, |§. 
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Figure 12: Predictions for rapidity distributions of several particles in p-p and central (b < 2 
fm) Au-Au collisions at RHIC from VNI, HSD and VNI+HSD models, taken from [||]. 



